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ABSTRACT 

Cooling flow models for the hot gas in elliptical galaxies predict that gas is cooling at a rate 
of 1 Mq yr~^, yet there is little evidence for this phenomenon beyond the X-ray waveband. If 
hot gas is cooling, it will pass through the 3 X 10^ K regime and radiate in the O VI AA1032,1038 
ultraviolet lines, which can be detected with the Far Ultraviolet Spectroscopic Explorer (FUSE) 
and here we report on FUSE observations of the X-ray bright early-type galaxies NGC 1404 and 
NGC 4636. In NGC 1404, the O VI doublet is not detected, implying a cooling rate < 0.3 Mq 
yr~'^, which is below the predicted values from the cooling flow model of 0.4-0.9 Mq yi~^. In 
NGC 4636, both O VI lines are clearly detected, indicating a cooling rate of 0.43 ±0.06 Mq yr~^, 
which falls within the range of values from the cooling flow prediction, 0.36-2.3 Mq yr~^and is 
closest to the model where the production of the cooled gas is distributed through the galaxy. 
The emission line widths, 44± 15 km s^^, are close to the Doppler broadening value (30 km s~^), 
indicating that the flow is quiescent rather than turbulent, and that the flow velocity is < 30 km 
s-\ 

Subject headings: galaxies: individual (NGC 1404, NGC 4636) — galaxies: ISM 



1. Introduction 

In the standard model for cooling flow ellip- 
ticals, the gas shed by stars during normal stel- 
lar evolution becomes thermalized, colliding with 
ejecta from other stars and with the ISM. This 
process heats the gas to lO^'^-lO^ K, and if there 
were no additional heating or cooling, the gas 
would be bound to the galaxy and have the same 
spatial distribution as the stars. However, ra- 
diative cooling drains energy from the gas most 
rapidly at small radii, causing a loss of buoyancy 
and a subsequent inflow of gas. The rate at which 
gas cools and flows inward (A/) is proportional 
to two observed quantities, the energy loss rate 
{Lx) divided by the thermal energy per gram 
likTx/tJimp), or M oc Lx/Tx, typically 0.03-3 
Mq yr~^ for ellipticals. 

Although cooling flow models do an adequate 
job in explaining the X-ray data, there is little 
conflrmation for the consequences of the models. 
In particular, the cooling flow will produce cooled 



gas, yet gas at < 10^ K is not seen in abundance 
in these systems (Roberts et al. 1991), suggesting 
that the gas is turned into stars. However, if stars 
are produced with a normal initial mass function, 
young blue stars would have been detected if the 
rate was above 0.01 Mq yr~i (O'Connell 1999). 
Therefore, either star formation in these systems 
does not produce many stars above 1 Mq, or the 
basic model is wrong. Most of the modifications 
to the basic model involve some sort of heating 
mechanism to balance the radiative cooling and 
prevent cool gas from being produced. The crit- 
ical issue is whether the gas is losing its thermal 
energy and a simple test is to search for the emis- 
sion lines that are produced as gas cools through 
intermediate temperatures. 

The most powerful diagnostic line for such a 
test is the O VI doublet, which is produced as 
the gas passes through the 3 X 10^ K range, at 
which point the gas has lost 90-95% of its orig- 
inal thermal energy, so it has effectively cooled. 
As gas cools through the 2-4x10^ K range, the 
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cooling is carried by a single ionization state of 
a single element, O VI, and its doublet lies at 
AA1032,1038, now accessible with the Far Ultra- 
violet Spectroscopic Explorer (FUSE). The rela- 
tionship between line luminosity and cooling rate 
(M) is insensitive to the metallicity of the gas or 
to whether the gas is out of equilibrium (Edgar 
and Chevalier 1986). Consequently, the luminos- 
ity of the AA1032,1038 lines is a direct measure of 
M. We have a program using FUSE (Moos ct al. 
2000) to observe cooling flow ellipticals and here 
we report on the first two observation, which are 
of the classic cooling flow galaxies NGC 4636 and 
NGC 1404. 

2. Observations 

Out of a more extensive program, NGC 1404 
and NGC 4636 are the first two galaxies to be ob- 
served, and they are X-ray bright, so their emis- 
sion is dominated by hot gas (e.g., Canizares, Fab- 
biano, and Trinchieri 1987). Furthermore, the 
cooling time of the X-ray emitting gas is far less 
than a Hubble time, making these classic cooling 
flow galaxies. There are important differences be- 
tween these two galaxies in their optical luminosity 
as well as their X-ray to optical luminosity ratio, 
leading to different values for the mass cooling rate 
M. The galaxy NGC 1404 is in the Fornax clus- 
ter and for a distance of 28.4 Mpc {Hq = 50 km 
s~^ Mpc~^; the conclusions are distant indepen- 
dent), logLs = 10.74, logLx = 41.27 and Tx = 
0.56 kcV (Faber et al. 1989; Brown and Bregman 
2000), while NGC 4636 is in the outskirts of the 
Virgo cluster with a distance of 26.7 Mpc and with 
properties log Lb = 10.96, logLx = 41.81 and Tx 
= 0.72 keV. 

The observations for NGC 1404 were obtained 
in 10 December 1999, with an exposure time of 
7.49 ksec, while the data for NGC 4636 were ob- 
tained in 23 May 2000 with an exposure time of 
6.46 ksec. In both cases, the large aperture was 
used, which is 30" square, leading to a velocity 
resolution of 100 km s~^; there is no effective spa- 
tial information perpendicular to the dispersion 
axis. The data streams did not contain periods of 
bursts or other detrimental events, so editing of 
the event file was not necessary. The data were 
processed with the pipeline program during Jan- 
uary 2001, including an optimal extraction algo- 



rithm, which uses a local background and employs 
information about the shape of the spectrum per- 
pendicular to direction of the spectrum. Four dif- 
ferent detector segments cover the spectral region 
of interest, two LiF segments and two SiC seg- 
ments, but the LiF la channel had the highest S/N 
and was used in most of the analysis. Addition of 
the other channels generally reduced the S/N and 
led to spurious features in the spectrum. The 30" 
aperture was centered on the optical centers of the 
galaxies and the resulting spectra around the re- 
gions of the O VI emission in NGC 4636 and NGC 
1404 are show in Figures 1 and 2. 

These spectra have features not associated with 
the elliptical galaxies, such as airglow lines and 
Galactic absorption lines. These lines of sight are 
out of the plane of the Milky Way, where the H I 
columns are 1.36 x 10^° cm-^ (NGC 1404) and 
1.63 X 102" (NGC 4636), and at those low 

values, the H2 absorption lines are not expected 
to be large. The strongest Galactic atomic line in 
the region is that of C II at 1036.34 A, followed 
by O I lines at 1026.47 A and 1039.23 A. A few 
other lines fall in this window, the ground-state 
doublet of O VI AA1031.95, 1037.63 and the ex- 
cited metastable line C 11* A1037.02, and although 
Galactic absorption by these lines is seen in a num- 
ber of other sightlines through the galaxy (e.g., 
Sembach et al. 2000), they are weaker. These two 
O VI lines are the ones predicted to be present 
in cooling fiows in emission, and for optically thin 
conditions, the O VI A1032 line will be twice as 
strong as the O VI A1038 fine. 

2.1. NGC 1404 

Of the four segments, only the LiFla had a 
count rate above the background rate and it pro- 
duces the highest S/N spectrum in the 1025-1050 
A range (Fig. 2). The spectrum in the range 1025- 
1050 A range shows a stellar continuum upon 
which terrestrial airglow lines plus a few Galac- 
tic absorption features are present (Fig. 2): the 
strong line of C II A1036.34 at zero velocity; an 
O I A1026.47 line adjacent to the Ly/3 airglow fea- 
ture; and Ly/3 absorption, although it is filled in by 
the Ly/3 airglow line plus two other airglow lines 
at 1027-1028 A. Given the modest intensity of the 
airglow line at 1025.5 A, the other airglow lines in 
this spectral region should be unimportant, and 
they are generally not visible (e.g., the O I A1039.2 
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line is absent). The strongest H2 lines in this and 
other speetral regions are not visible, so the Galac- 
tic N{H2) < 10^^ cm~^. There are no interstellar 
atomic or molecular absorption lines at the red- 
shift of NGC 1404. 

The stellar continuum has one dominant fea- 
ture in Figure 2, Ly/3 absorption, plus a number 
of minor absorption lines. In an effort to under- 
stand which features are stellar and which might 
be interstellar, we can compare it to an approriatc 
stellar model. As discussed by O'Connell (1999), 
the continuum in this region is produced by low 
mass stars that cither populate the extreme hori- 
zontal branch or have evolved away from it. Mod- 
els that include these various stellar contributions 
were calculated by Brown et al. (1997) for com- 
parison with HUT data, with atmospheric metal- 
licities of O.IZq leading to the best fits. These 
low atmospheric abundances probably reflect dif- 
fusion processes since the mean metallicities of the 
stars in these systems are near-solar. The model 
of Brown et al. (1997) is at 3 A resolution, which 
is too low for comparison with our data, so instead 
we use the average of two stars of different grav- 
ity {9 = 4,5.5), as these are representative of the 
range of gravities for the stars that contribute to 
the flux in this wavelength region (Dorman, Rood, 
and O'Connell 1993); we used a surface metallic- 
ity of O.IZq and T = 26,000 K. Also, this model 
was convolved with the velocity dispersion of the 
galaxy (Fig. 2). 

The expected O VI emission lines from NGC 
1404 would be shifted to AA1038. 6, 1044.3 A(1947 
km s~^), and relative to the continuum, there is no 
evidence of emission from cither line. To estimate 
the upper limit to the line strength, we assume 
that the line width is characteristic of the veloc- 
ity dispersion of the system, AA = 1.0 A (300 
km s~^). Conservative 3(t upper limits for the 
lines are ^'(A1032) < 3 x lO'^^ erg cm^^ s^^, and 
F(A1038) < 1.4 X 10-1^ erg cm^^ s^^. The ex- 
tinction in this region is very low, estimated to be 
Ay = 0.00-0.05, so the average would lead to a 
extinction correction of a factor of 1.08 at 1035 A. 

2.2. NGC 4636 

The spectrum of NGC 4636 contains two Galac- 
tic features seen in the spectrum of NGC 1404, the 
C II A1036.34 and the O I A1039.23 feature (Fig. 
1), although they are shifted by 0.24 A and 0.30 



A respectively from their rest wavelengths (about 
-|-70 km s~^). These Galactic absorption lines are 
expected to occur near km s~^, yet there is no 
line at this location with the instrumental width. 
Also there is no known Galactic H I feature of 
any importance at -|-70 km s~^, which would need 
to dominate the neutral column in this direction 
to cause the observed results. Instead, it appears 
that there is a spectral shift of about 0.3 A, and a 
similar shift is seen in the Ly/3 airglow line, which 
is redshifted by 0.40 A. Consequently, in our anal- 
ysis, we apply a shift to the data of -0.3 A (—87 
km s"'^). In addition to the atomic lines, there 
is modest absorption by H2, although only the 
strongest lines are detected, such as the 1049.7 
A line. Based on the strength of this and other 
Galactic H2 (e.g., in the 1062.5-1064 A region), 
we estimate that N{H2) » 10^^-^ cm"^, typical of 
lines of sight with this 7V(H I) value. 

Aside from the absorption lines and airglow 
lines, the two strongest features are emission lines 
that correspond to the wavelengths at which the 
O VI doublet would occur near the redshift of 
NGC 4636; the line centers are 1018 ± 20 km s"^ 
(A1032 line) and 1061 ± 20 km s^^ (A1038 line) 
compared to the galaxy redshift of 938 ± 4 km s~^ 
(from NED). The separation of the line peaks is, 
to within the errors, consistent with the separa- 
tion of the two lines in the O VI, but the mean 
velocity of the lines is 100 km s~^ greater than 
that expected. This difference may result from 
wavelength calibration errors since in a previous 
version of the pipeline processing, the wavelength 
of the O VI was consistent with optical redshift for 
NGC 4636. 

The two lines are single-peaked and fairly nar- 
row. The FWHM of the Unes are about 0.15±0.05 
A (44 km s^^), with the error depending upon 
where one defines the continuum. The instrumen- 
tal resolution depends upon the degree to which 
the aperture is illuminated. The airglow lines, 
which fill the aperture uniformly, have a FWHM 
of 0.34 A (100 km s^^), while a point source would 
have a FWHM of about 0.05 A (15 km s^^). 
Therefore, the emission lines from NGC 4636 must 
not fill the aperture uniformly, although we can- 
not determine whether the line width is due to an 
instrumental limitation or to Doppler broadening. 
The line fiuxes are J^(A1032) = 3.0 ± 0.6 x lO^^'"^ 
erg cm-2 s"!, and F(A1038) = 2.4 ± 0.6 x 10"^^ 
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erg cm"^ s~^, with the primary error in the Hne 
fluxes resulting from the uncertain placement of 
the continuum. The line ratio, 1.25 ± 0.32, is 
2.3(7 from the optically thin value of 2.0, but 
given the uncertainty in the model stellar contin- 
uum, we do not consider this difference particu- 
lary significant. The extinction along this sight- 
line is Ay = 0.04-0.09, depending upon the red- 
dening model used (Burstein and Heiles 1982, or 
Schlcgcl, Finkbcincr, and Davis 1998), so we use 
the average of the two values, along with a value 
of ^(1035 A)/Av = 4.0 (CardelU et al. 1989), 
which leads to an cxctinction-corrccted value of 
f (A1032 + A1038) = 6.8 ± 0.9 x lO'^^ erg cm'^ 

3. Interpretation and Discussion 

The conversion between the O VI line luminos- 
ity and cooling rate depends only weakly on the 
details of the physical situation, provided that the 
initial temperature is significantly higher than the 
temperature at which the ion is most abundant. 
This situation occurs in these galaxies, where the 
ambient temperature of the hot X-ray emitting gas 
is 6-9x10^ K and the ionization fraction in O VI 
is negligible. The cooling of the gas will proceed 
isobarically if the cooling time is long compared to 
the sound-crossing time of the region, eventually 
becoming isochoric as the gas cools and the rela- 
tive size of the timescales is reversed. If gaseous 
cooling occurs in regions smaller than 1 kpc, the 
transition from the isobaric to isochoric case oc- 
curs at T < 4 X 10^ K (Edgar and Chevalier 1986), 
and the cooling may be entirely isobaric as it ap- 
plies to O VI. Unfortunately, we do not at this time 
have a good contraint on the size of the cooling re- 
gion. For the pure isobaric case, Edgar and Cheva- 
lier (1986) show that L(A1032) = 1.1 x lO^^M 
erg (M is in Mq yr~^) and for the pure iso- 
choric case, L(A1032) = 0.7 x lO^^M erg s'^ (the 
A1038 line is half of the strength of the A1032 
line). In this paper, we will use a conversion that is 
halfway between the isochoric and isobaric cases, 
i(A1032) = 0.9 X lO^^M erg s-\ which introduces 
only a 20% uncertainty in the model; the resulting 
values of M are given in Table 1. 

The critical test is to compare these values of M 
with those derived from X-ray measurements. The 
cooling rate M is the rate at which the thermal en- 



ergy of the system is being drained away by radia- 
tive losses, so it is the ratio of the net cooling rate 
divided by the specific thermal energy content, or 
M — Lnet/E, where the specific thermal energy is 
E = 3kT/ finip, jjL is the mean molecular weight for 
the hot gas (/x = 0.63) and rup is the usual proton 
mass. The net energy loss rate Lnet equals the 
radiative loss rate (the bolometric X-ray luminos- 
ity) minus heating sources such as supernovae or 
gravitational compression as the gas falls inward. 
The amount of gravitational compressional heat- 
ing depends upon whether the gas cools close to 
where it entered the flow or whether it flows into 
the central region (r < Tcore) before losing most of 
its thermal energy. The gas will become cool only 
in the central region unless thermal instabilities 
can grow throughout the flow, but the best cal- 
culations on this matter indicate that most linear 
perturbations will not grow rapidly enough to be 
effective (Balbus 1988). However, models without 
distributed mass drop-out lead to an X-ray sur- 
face brightness profile that is too sharply peaked 
in the center, inconsistent with observations. Con- 
sequently, most models include mass drop-out as 
a function of radius and the usual formulation is 
that the rate of mass drop-out is inversely propor- 
tional to the cooling time, ov p = qp/tc, where tc is 
the isobaric cooling time (Sarazin and Ashe 1989). 

Without mass drop-out, q = 0, and a moder- 
ately efficient rate of drop-out would be given by 
q = 1, so we consider both cases. We use the 
model from Sarazin and Ashe (1989), where Lb = 
1 X 10^^ Lq, and expressing Lx in units of 10^^ erg 
s~i and Tx in keV, we obtain M = 0.26LxT^^ 
Mq yr-1 for g = 0, and M = OAOLxT'^ Mq 
yr~^ for g = 1 (this is for the absorption-corrected 
X-ray luminosity in the 0.5-2.0 keV band). In the 
q = 1 case, the mass drop-out occurs throughout 
the galaxy and since our 30" square aperture does 
not encompass the entire galaxy, we must deter- 
mine M within the aperture. For our two galaxies, 
the square aperture encloses the emission from the 
galaxy within about 2.2 kpc (NGC 4636) and 2.4 
kpc (NGC 1404) from the center, and in the q = 1 
model, this would enclose about half of the total 
M. After applying this factor to the q = 1 case, 
the resulting M is similar to the q = case, with 
M = 0.20LxTx^ Mq yr"^. We use the values of 
Lx and Tx given in Brown and Bregman (2000) 
for these galaxies to calculate values of M for the 
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q — case and for the q = I case for the coohng 
only within the aperture (denoted Mo, Mi and 
given in Table 1). 

An alternative method of determining the ex- 
pected value of M is to use the X-ray observations 
directly by noting that p = qp/tc oc p^/A(Tx), but 
since the cooling function A{Tx) changes far less 
than p in these galaxies, a good approximation is 
p oc p"^, and the amount of cooling gas in the aper- 
ture is the product of the line integral of p through 
the galaxy, or M ex ^ / p^dl, where A is the area of 
the aperture. The quantity / p^dl is the emission 
measure, to which the X-ray surface brightness is 
proportional, so the fractional amount of M within 
the aperture is just MA/Mrot = Lx,A/Lx,Tot- 
We have calculated the quantity Lx,a/ Lx.Tot di- 
rectly from the ROSAT PSPC data for these sys- 
tems (using the backgrounds discussed by Brown 
and Brcgman 2000) and obtain the fractional cool- 
ing within the FUSE aperture to be 0.083 (NGC 
4636) and 0.25 (NGC 1404). This leads to values 
of M within the apertures of 0.36 M© yr^^ (NGC 
4636) and 0.40 M© yr"! (NGC 1404), which is 
the product of Lx. a/ Lx.Tot and the value of M 
given above for the q ~ \ case (given as M2 in 
Table 1). For NGC 1404, this value of the cooling 
rate in the aperture is within a factor of two of 
that derived entirely from the Sarazin and Ashe 
model. However, in the case of NGC 4636, the 
two values differ by up to a factor of six. Prob- 
ably, this difference arises because NGC 4636 is 
more extended than the average galaxy of simi- 
lar optical luminosity. Whereas NGC 4636 and 
NGC 1404 are similar in both optical luminosity 
and distance, NGC 4636 has a value for the half- 
light radius, r^, that is nearly four times greater 
than for NGC 1404 (101" compared to 27"). Con- 
sequently, the core radius used in the model by 
Sarazin and Ashe (1989) may be too small in the 
case of NGC 4636 and this would lead to an over- 
estimation of the fractional amount of M within 
the aperture for the q = 1 case; the values of M 
are not affected for the q — Q case. The different 
values for M from the X-ray data may be regarded 
as the expected range in the prediction, given the 
uncertainties in the models. 



4. Conclusions 

For NGC 4636, the value of M from the O VI 
measurement lies within the predicted range of M 
from the X-ray data, lending strong support to 
the cooling flow model. Furthermore, the O VI 
result indicates a cooling rate significantly below 
that predicted from the q — Q model, a discrep- 
ancy that is removed if wc assume that the cooling 
is distributed throughout the galaxy, which also 
was predicted from the application of the cooling 
flow models to the X-ray surface brightness dis- 
tribution. Since M(0 VI) « M2, the mass drop- 
out is consistent with that obtained by assuming 
that M(< r) oc Lx{< r), although more detailed 
models providing determinations of M should be 
calculated for NGC 4636, such as along the lines 
of Bringhenti and Mathews (1996, 1999), who in- 
clude a more accurate galaxy model plus impor- 
tant effects of accretion within the galaxy group. 

In NGC 4636, the O VI line width of 44 ± 15 
km is similar to the thermal Dopplcr width 
of 29 km s~^ and is less than either the sound 
speed in the gas, 77 km s~^, or the velocities of 
the 10^ K optical emission line gas of 200 km s"-"^ 
(Caon et al. 2000). The narrow line width has 
several implications, the first being that the radial 
velocities within the cooling fiow arc less than 30 
km s^^ (3 a upper limit). This is consistent with 
the cooling flow calculations, which predict typical 
flow velocities of ^ 10 km s^^ (e.g., Sarazin and 
Ashe 1989). Furthermore, the environment of the 
O VI gas must be quiescent, also consistent with 
models. 

In contrast to NGC 4636, for NGC 1404 the 
predictions of the cooling flow model are not con- 
firmed in the O VI data, where our upper limit 
is below the lowest X-ray prediction for M. The 
most likely resolution of this discrepancy is that 
the net radiative loss rate is not given by the X- 
ray luminosity, which will occur if there is an addi- 
tional source of heating, such as from supernovae 
or a central AGN. There is no evidence for a cen- 
tral radio source, as there is only an upper limit 
to the radio continuum emission of 0.7 mJy (6 cm; 
Sadler, Jenkins, and Kotanyi 1989). 

In the future, we will observe additional galax- 
ies to determine whether the cooling flow model 
is generally applicable. Also, we will obtain ad- 
ditional observations of NGC 4636 in off-central 
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locations, in order to determine the radial distri- 
bution of the cooling material, which is vital to 
modelers and to the determination of the location 
of the deposition of cooled gas that may produce 
new stars. 
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Table 1 

Cooling Rates from X-Ray and O VI Observations.'^ 



Galaxy Name 


M(0 VI) 


Mo 


Ml 


M2 


NGC 1404 


< 0.3 


0.86 


0.66 


0.40 


NGC 4636 


0.43 ± 0.06 


2.3 


1.8 


0.36 



''Units of M are Mq yr~^ 




1030 1040 1050 

Wavelenth (A) 



Fig. 1. — The FUSE spectrum of NGC 4636 in the 1020-1050 A region shows airglow hnes (notably Ly/3; 
denoted by +), several Galactic absorption lines, the strongest being C II (1036.34 A) and O I (1026.47 A, 
1039.23 A), and possibly O VI (1032 A), as well as the redshifted O VI emission lines of from NGC 4636 
(in italics). The dashed line is the expected stellar continuum from NGC 4636. The spectrum has been 
smoothed to a velocity resolution of 100 km s~^ for display purposes, although from the full resolution data, 
the O VI emission lines have a FWHM of 44 km s~^, close to the thermal width of 30 km s~^. 
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1030 1040 1050 

Wavelenth (A) 



Fig. 2. — The FUSE spectrum of NGC 1404, similar to that in Figure 1, except there is no apparent detection 
of the redshifted O VI emission lines from NGC 1404. 
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